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Moduli dark matter and the search for its decay line using Suzaku X-ray telescope 
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Light scalar fields called moduli arise from a variety of different models involving supersymmetry 
and/or string theory; thus their existence is a generic prediction of leading theories for physics 
beyond the standard model. They also present a formidable, long-standing problem for cosmology. 
We argue that an anthropic solution to the moduli problem exists in the case of small moduli masses, 
and that it automatically leads to dark matter in the form of moduli. The recent discovery of the 
125 GeV Higgs boson implies a lower bound on the moduli mass of about a keV. This form of 
dark matter is consistent with the observed properties of structure formation, and it is amenable 
to detection with the help of X-ray telescopes. We present the results of a search for such dark 
matter particles using spectra extracted from the first deep X-ray observations of the Draco and 
Ursa Minor dwarf spheroidal galaxies, which are dark matter dominated systems with extreme mass- 
to-light ratios and low intrinsic backgrounds. No emission line is positively detected, and we set 
new constraints on the relevant new physics. 



I. THE NATURE AND ORIGIN OF MODULI 

Superstring theories, widely considered as candidates 
for a unified theory of all interactions, generically pre- 
dict the existence of light scalar fields associated with 
the breaking of scale invariance and the size and shape 
of the string compactification volume. One example is 
the dilaton, a scalar field associated with the scale trans- 
formations whose vacuum expectation value (VEV) de- 
termines the values of various couplings in the low-energy 
effective field theory. String theory compactification radii 
also appear in the low-energy theory as scalar fields with 
some very small masses. In a number of models, the 
string scale M s is of the order of the reduced Planck 
mass, Mq = M P /V87T = 2.4 x 10 18 GeV, but, in models 
with a large compactification volume, the scale of M s is 
suppressed by the compactification volume and can be 
much lower [l| . Although one might expect that interac- 
tions of moduli with other fields should be suppressed by 
a factor that is proportional to some power of Mq, de- 
tailed calculations can produce answers that are different 
from naive dimensional analyses Q. 

Supersymmetry is a generic prediction of string theory, 
but, even aside from string theory, it is a well-motivated 
concept in its own right. An appealing generaliza- 
tion of space-time symmetries involving non-commuting 
fermionic degrees of freedom in the form of a graded alge- 
bra leads to supersymmetry and supersymmetric exten- 
sions of the standard model. As long as supersymmetry 
is unbroken, its potential has numerous flat directions 
with infinitely many degenerate classical vacua. The de- 
grees of freedom parameterizing these flat directions are 
massless in the limit of unbroken supersymmetry, but 
they acquire a mass from the breaking of supersymme- 
try. The mass thus acquired is well below the scale at 
which supersymmetry breaking occurs. For example, for 
-F-type breaking, the breaking of supersymmetry by the 



non-zero VEV of \F\ can be communicated to the rest of 
the fields by gravitational interactions suppressed by the 
(reduced) Planck mass Mq. This gives moduli masses of 
the order of m ~ \F\ /Mq, well below the scale of super- 
symmetry breaking. The coupling of the moduli to the 
rest of the fields is also suppressed by the higher scale. 



II. THE MODULI MASSES 

Massless in the limit of exact supersymmetry, the mod- 
uli acquire some masses from supersymmetry breaking. 
The range of moduli masses is model-dependent, and a 
lighter scalar with mass well below a keV would not be 
a viable dark matter candidate. However, the recent dis- 
covery of a 125 GeV Higgs boson [3, [j] points to the 
moduli masses above 1 keV, which further strengthens 
the motivation for investigating moduli as a viable dark- 
matter candidate. 

In the following section we directly address the moduli 
problem, and our approach is based on a low-scale super- 
symmetry breaking. Therefore, we consider the gauge- 
mediated supersymmetry breaking scenario. The break- 
ing of supersymmetry occurs via a non-zero vacuum ex- 
pectation value (VEV) of the F-component (Fs) ^ of 
some chiral superfield S, whose scalar component also 
has a non-zero VEV (S). Supersymmetry breaking is 
then communicated to the visible sector by messengers 
v&i. The messengers couple to S with couplings Ay, 



W = A,,.S>lvl/, 



(1) 



In what follows we suppress the indices as we discuss the 
generic features of a class of models. The mass-squared 
matrix of the scalar messengers takes the form 
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The stability condition requires that this matrix have no 
negative eigenvalues, which implies 



M, 



2 _ 



|A(5)| 2 > \X(F S )\. 



(3) 



In the visible sector, the squarks acquire masses through 
the gauge interactions involving the messengers in the 
loops. Up to some group theoretical factors of order one, 
these masses are 



Q3 MFs) 
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(4) 



where a 3 is the gauge coupling constant. The recent 
discovery of the Higgs boson with mass 125 GeV sig- 
nificantly tightens the allowed range of parameters for 
the squark masses. According to the analysis of 0, E3], 
a mass of 125 GeV requires that the squarks be heavier 
than 10 TeV. The constraint 



m sq > 10 TeV, 



(5) 



combined with Eq. the stability constraint Eq. 
and the requirement that A < 1 implies the following 
lower limit on the supersymmetry breaking scale = 
\Fs \ + (other contributions) > \Fs\- 



\F\ > \Fs\ > 



10 TeV 



^10 6 GeV)' 



(6) 



We note that \F\ = \Fs\ in the case of direct media- 
tion 0. 

This lower bound on the supersymmetry breaking scale 
implies the following lower bound on the mass of the 
moduli: 



J£l 
M G 



> 1 kcV. 



(7) 



Some additional constraints arise from the require- 
ment of not over-producing dark matter in the form of 
gravitinos, which have a similar mass. For masses in 
the several keV range the gravitinos come into thermal 
equilibrium and reach thermal abundances at temper- 
atures above TeV, and the total amount of dark mat- 
ter would be unacceptably large if no dilution occurred. 
However, gauge-mediated supersymmetry breaking mod- 
els contain candidate particles whose out-of-equilibrium 
decay can produce entropy and dilute the population 
of dark matter to an acceptable level || @. The re- 
quired dilution factor A is proportional to the mass, 
A ~ 100(TO3/2/10keV) Q. For masses above 1 McV, 
the requisite dilution of A ~ 10 4 can be difficult to rec- 
oncile with leptogencsis because the baryon asymmetry 
is diluted by the same factor. For masses below (but 
not too far below) 1 MeV the dilution can still be ef- 
ficient, and the gravitino abundance well below the ob- 
served dark matter abundance. As we discuss below, 
the moduli generally tend to be over-produced, but an- 
thropic selection considerations point to abundances in 
the range of observational estimates of the dark matter 



density. Since the same anthropic reasoning would not 
work for the gravitinos, one can assume that the mod- 
uli have a greater abundance than the gravitinos. It is 
also possible that the contributions of the moduli and 
the gravitinos are comparable. The intriguing possibility 
that dark matter is composed of two separate compo- 
nents with different free-streaming properties, which can 
have some observable effects on the small scales, was dis- 
cussed in detail in connection with sterile neutri- 
nos produced by two distinct mechanisms [TFj, [JJ|-[l|| . In 
our case, there is an analogous possibility that dark mat- 
ter is comprised of two components, namely, the moduli 
and the gravitinos. 

Taking into account the Higgs boson mass and the 
gravitino constraints, the moduli masses of interest lie 
in the range ~ (1 — 10 3 ) keV. 



III. THE MODULI PROBLEM 

As discussed above, a number of independent argu- 
ments based on well-motivated theories converge on the 
prediction of light scalar fields coupled very weakly to 
the remaining fields in the theory. The relevant masses 
range from a few eV to well above the electroweak scale. 
The lighter moduli can decay through interactions sup- 
pressed by an effective high energy scale A c ff into photons 
and neutrinos. Typically, the dominant decay channel is 
(j) — > 77 through a coupling 

where we have defined a parameter b, 



b = 



A of f : 



(9) 



choosing b and A c fj so as to parameterize the coupling 
in a model-independent way. A s and the correspond- 
ing parameter b have specific meanings in the context of 
specific models. For example, for string moduli, A e ff can 
be related to the string compactification scale and struc- 
ture. Explicit calculations yield values of b ranging from 
b = \/2 ~ 1 [rjj to much larger values . 
The rate of moduli radiative decay is 



r — ^2 



(10) 



The generic prediction of light scalar fields gives rise to 
a very serious cosmological problem, namely the moduli 
problem. There is strong and growing evidence that the 
universe has undergone a period of inflation that resulted 
in flatness and homogeneity of the observed universe, and 
that generated the observed spectrum of density pertur- 
bations. If mass 77i0 is smaller than the Hubble constant 
H during inflation, then the vacuum expectation value 
of the light scalar field during inflation can be very large 



3 



for the following reason. One expects the energy den- 
sity in each scalar field to be of the order of H 4 , that is 
H 4 , during inflation. Hence, the VEV (<fi), on 



m 



2^,2 



average (i.e., averaged over superhorizon scales), takes a 
very large value. Furthermore, the value of the ground- 
state VEV is not well defined. During inflation, the min- 
imum of the potential for the <j) field can differ from its 
zero-temperature value due to the terms associated with 
supersymmctry breaking in de Sitter space, i.e., due to 
expansion of the universe. These terms, coming from the 
Kahler potential, act as additional mass terms cH 2 <fr 2 , 
where c is a constant of order 1 that can be positive 
or negative [H, [lj| . These contributions overwhelm the 
small moduli masses and can displace the VEV of the 
field by a large amount. When inflation is finished, and 
reheating commences, the VEVs of the moduli fields are 
stuck at some values far from the minima of their respec- 
tive potentials. Eventually the Hubble constant H ap- 
proaches the mass m^; and, only then can the field start 
oscillating. However, the minimum of the effective poten- 
tial, which includes the expansion-induced contributions, 
can still differ from the zero temperature minimunQ. 

A scalar field with some vanishingly small couplings to 
the lighter particles oscillates until the time 



= r^ 77 = 7.6 x 10 



:S2 



1\ 2 /lkeV x 3 



(11) 



During this time, the field's cosmological behavior is 
identical to that of non-relativistic matter, with energy 
density declining with scale factor a as 1/a 3 . When the 
radiation density that scales as 1/a 4 drops below the en- 
ergy density of moduli, the universe enters a prolonged 
period of matter dominated expansion, at the end of 
which moduli decay. This causes a variety of cosmolog- 
ical problems, depending on the moduli mass and decay 
width. 

If moduli decay before big-bang nucleosynthesis, their 
decay products thermalize and they do not cause a con- 
flict with observations, assuming that the baryon asym- 
metry of the universe is not diluted away. It is possible 
to generate the baryon asymmetry through Affleck-Dine 
baryogenesis j20l - |22| so as to compensate for dilution. 
Masses above 10 2 TeV may be considered "safe" in 
this sense. If the moduli decay during big-bang nucle- 
osynthesis, the entropy release may still be consistent 
with successful nucleosynthesis, and the additional en- 
tropy released may account for the small observed de- 
viation of the effective number of light species from the 
standard model prediction [2^, [24|. This case requires 
further investigation. 



1 This stymies some of the attempted solutions: even if the mod- 
uli oscillations can be damped at some point in the expanding 
universe for some specific value of the Hubble parameter H, they 
are bound to resume when H decreases further because the min- 
imum of the effective potential depends on H [TJJ . 



Moduli decay after nucleosynthesis, but before recom- 
bination, unacceptably dilutes the baryon density. Decay 
after recombination also distorts the cosmic diffuse mi- 
crowave, and contributes to the X- and gamma-ray, back- 
grounds in violation of observational constraints for some 
combinations of abundance and mass [25U27I ] . although 
cosmologically significant densities are not excluded for 
Tn^ < 100 keV Um. 

Finally, although decay times longer than the age of 
the universe are consistent with the moduli being the 
dark matter if they can be produced with the correct 
abundance, naive estimates for moduli dark matter abun- 
dance in this case predict a much greater value than the 
observed dark matter density. Indeed, let us estimate 
the energy density of the universe when the modulus <p 
starts oscillating, that is, when H ~ m^, or when the 
temperature of the universe is 



T ~ {QO/iPgjWy/MGmi. 



(12) 



The density to entropy ratio is 

s ~ (27r 2 /45) 5 *T| ~ ^ VkcVJ \M G 

(13) 

Comparing this ratio, corrected for the entropy produc- 
tion, with today's value, 



Pdm 



0.2^ 



3 x 10- 10 GeV, 



(14) 



clearly reveals a drastic discrepancy. 

With the help of thermal inflation [29l[30|. or by means 
of coupling the moduli fields to the inflaton [HI, HH the 
moduli problem can be ameliorated in some range of 
masses [2a - l27j . However, neither of these approaches 
provides a complete solution. 



IV. ANTHROPIC SELECTION 

As long as the moduli are long-lived, anthropic selec- 
tion may be invoked to explain why moduli dark matter 
in the observed part of the universe is present in the right 
amount, circumventing the problem discussed in the pre- 
vious section of the energy density in the oscillating mod- 
uli field taking on an excessively large value. The moduli 
problem is eliminated if the initial value </>o of the mod- 
ulus field is, for whatever reason, sufficiently small that 
the energy density in the moduli does not significantly 
exceed the baryonic energy density. This may be com- 
pared with the requirement, realized in the observable 
universe, for the dark-to-baryonic matter ratio to be ~ 4 
in order for conditions amenable to the the formation 
of stars and planets, and, ultimately, for the existence 
of life, to be attained [331 ] . Therefore, one may simul- 
taneously attempt to explain away the moduli problem 
and explain the abundance of dark matter based on the 
anthropic selection: life can only emerge in those parts 
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of the universe where the moduli energy density is suf- 
ficiently low that the dark-to-baryonic matter ratio has 
the acceptable value. 

The reasoning here is analogous to that applied to the 
axion (34[. Indeed, the abundance of string moduli de- 
pends on the initial value <j>o (Eq. [T3]) that is determined 
by the VEV of the field during inflation. On average, 
this VEV tends to be very large, of the order of Mq. 
However the large VEV may generate too much matter 
that is already dominant at a temperature T ~ 10 5 GeV. 
Inflation still assures that fi = 1, so that the universe 
remains flat; however, in this case structure formation 
proceeds very early leading to a universe dominated by 
black holes [33|, |35[. If matter comes to dominate the 
energy density of the universe too early, the density per- 
turbations can grow and become non-linear before recom- 
bination. In this case, baryons and radiation are trapped 
together inside the collapsing halos, and the baryon cou- 
pling to photons maintains the Jeans mass at a constant 
value as the collapse proceeds. As a result, the bary- 
onic matter cannot fragment, and large amounts of cou- 
pled baryon-radiation fluid are dragged into the poten- 
tial wells created by clumps of dark matter. As discussed 
by [H| , the end result of such a structure formation pro- 
cess is a universe with supermassive black holes, pho- 
tons, and neutrinos, but without stars and planets. Such 
a universe is not amenable to life. For the example in 
question, the matter-radiation equality would occur at 
temperature 




Assuming that the baryon asymmetry escapes dilution, 
baryons are coupled to radiation until T rec ~ 1 eV. Lin- 
ear growth of density perturbations may commence at 
rcdshift z ~ 10 13 . Prior to recombination, even the over- 
densities as small as 10~ 9 should have gone non-linear. 
In the observed universe, primordial perturbations of the 
order of Sp/p ~ 10~ 5 are consistent with the data. Per- 
turbations as small as a one-millionth this level would 
have gone non-linear long before radiation decoupling. 
According to [35|, very small values of Sp/p < 10~ 8 are 
disfavored by some anthropic arguments as well. We con- 
clude that intelligent life is impossible in a universe where 
the dark matter density exceeds the baryon density by 
many orders of magnitude, and we refer the reader to 
[33l HH for a more detailed discussion. 

While the "most common" values of the initial field 
<p produce universes without stars or planets, inflation 
does allow for some relatively small part of the volume 
to have <j>Q <C Mq. The probability of this happening at 
any given point in de Sitter space is negligible, but the 
conditional probability of a small VEV, under the con- 
dition that life can evolve and produce an observer, is 
actually of the order of 1 because life is highly improba- 
ble in the large- VEV regions. Hence, anthropic selection 
would favor 0o ~ 10 _8 Mg, near the boundary of the 
parameter space acceptable for the existence of life. 



This motivates the search for dark matter in the form 
of string moduli with masses below 1 MeV and for decay 
widths corresponding to b = 1 — I0 3 . Cosmological data 
limits the lifetime of dark matter to T0_> 77 > 10 24 s in 
the case of radiative dec ays j36l . l37t | (and r 0—s. 77 > 3 x 
10 18 s for general decays [33]). Based on Eq. (fTTjl. this 
corresponds to sub-MeV masses. In particular, one can 
search for the dark matter in the form of moduli using 
X-ray telescopes. 



V. SEARCH FOR MODULI DARK MATTER IN 
X-RAYS 

The decay of moduli into two photons presents an op- 
portunity for a possible discovery using X-ray telescopes. 
Since this is a two-body decay, the signal is a very narrow 
line, whose width is determined by the Doppler broaden- 
ing due to the motion of the dark matter particles. The 
energy of a photon is (1/2) of the particle mass. If such a 
line is detected by the X-ray telescopes from compact as- 
trophysical objects with substantial dark mater content 
and/or from diffuse radiation due to the galactic and cos- 
mological distributions of dark matter, the wavelength of 
the line will provide a measurement of the dark-matter 
particle mass. The identity of the dark matter will still 
require some additional data, because moduli are not the 
only dark matter particles capable of producing an X-ray 
line. For example, sterile neutrinos, which can also be 
dark matter ^3Hl6l . |38| , are also expected to produce a 
line from their decay. Therefore, a discovery of an X-ray 
line from sources correlated with the dark matter distri- 
bution would reveal the mass of the dark matter particle, 
but not necessarily its identity. 



A. X-ray data used for this study 

As a proof of concept, we derive limits on the moduli 
decay rate in the soft X-ray energy band (~ 0.7—10 keV), 
based on observations of the Ursa Minor and Draco dwarf 
spheroidal galaxies conducted using the Suzaku X-ray 
Observatory. These data were obtained as part of a 
search for dark matter in the form of sterile neutrinos [39r - 
l4l| . The advantage of using Suzaku data, as detailed 
in [39l ] . lies in the low and stable internal background 
of the Suzaku XIS (X-ray Imaging Spectrometer) CCD 
detectors that, in combination with its sharp (by X-ray 
CCD standards) energy resolution, provides a robust and 
relatively sensitive capability for detecting faint spectral 
features over a broad (soft X-ray) bandpass. The moti- 
vation for targeting dwarf spheroidal galaxies is based on 
their proximity, high dark matter density, and absence 
of competing X-ray sources [3^, Ho[. The Ursa Minor 
and Draco systems, that have among the highest known 
dark matter surface densities (42| , remain the only dwarf 
spheroidals observed by Suzaku to date. 

Utilizing these data, we derive limits on T0_}. 77 , and 
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Galaxy Detector Time Total 


Source 


UMinor XIS03 


138.1 


10264 


5745 


UMinor XIS1 


69.05 


8949 


4380 


Draco XIS03 


124.8 


10361 


6396 


Draco XIS1 


62.4 


9173 


4636 



TABLE I. Exposures Times (in ks) and 0.6-7 keV Counts 



hence on the moduli b parameter, for in the 1.5 
20 keV mass range, as we now discuss in detail. 



B. Data reduction, spectral analysis, and line flux 
limits 

Spectral analysis of the Suzaku XIS data for Ursa Mi- 
nor and Draco were previously presented in [3^ and [!(| ■ 
Given subsequent significant enhancements in calibration 
data, and in reduction and analysis tools, we reprocess 
and reanalyze the data as follows, and also consider the 
datasets jointly for the first time. 

Our reprocessing proceeds along standard linct0 and 
follows that presented in j3^|, and more recently in [43| . 
utilizing the Suzaku aepipeline (version 1.0.1) ftool task. 
Details may be found in those papers (and references 
therein). Spectra are extracted from the entire fields- 
of-view for the back-illuminated (BI: XIS1) and two op- 
erational front-illuminated (FI: XIS0 and XIS3) chips, 
except for 4'-radius circular regions around the bright- 
est point source in each XIS field of view, and 2'-radius 
circular regions around two fainter sources in the Draco 
field. Spectral redistribution matrix files (rmf) are gen- 
erated using xisrmf gen version 2011-07-02, effective area 
function files (arf) using xissimarf gen version 2010- 
11-05. Spectra from the two FI detectors (XIS03 = 
XIS0 + XIS3) are co-added and a weighted XIS03 re- 
sponse function calculated from their respective rmf and 
arf files. Finally, non- X-ray particle background (NXB) 
spectra are extracted from the appropriately selected and 
weighted night earth data using xisnxbgen version 2010- 
08-22. All spectral fitting is conducted using Xspec ver- 
sion 12.70 Final good exposure times, and total (includ- 
ing NXB) and "source" (NXB-subtracted) counts are dis- 
played in Table 1. 

We adopt two approaches for establishing baseline 
models via spectral fitting. In the first the NXB is sub- 
tracted, spectra in the 0.6 — 7 keV bandpass are grouped 
into bins with a minimum of 15 cts, and best-fit mod- 
els are found by minimizing the x 2 statistic. The mod- 
els represent the astrophysical background, and are com- 
posed of Cosmic X-ray Background (CXB) and Galactic 
X-ray Background (GXB) components. The former is 



characterized as a power law, the latter as a soft thermal 
plasma. For Draco a two-temperature GXB model was 
required. Fits are performed jointly for Ursa Minor and 
Draco (but separately for XIS03 and XIS1). However, 
only the CXB slope is assumed identical in the two fields; 
GXB temperatures and all normalizations (expressed in 
units per solid angle) are untied. Reduced \ 2 m the 
best fit models are 0.97 (1038 degrees of freedom) for the 
XIS03, and 1.02 (918 degrees of freedom) for the XIS1 
fits. These spectra, best-fit models, and residuals are 
shown in Figure 1. 
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FIG. 1. 0.6 - 7 keV NXB-subtracted Suzaku XIS03 (left panel, 
(a)) and XIS1 (right panel, (b)) spectra of the Ursa Minor (black) 
and Draco (red) Dwarf Spheroidal Galaxies. The best joint-fit as- 
trophysical background baseline models and residuals (data-minus- 
model) are also shown. 



In the second approach, we consider the unsubtracted 
(total), unbinned spectra divided into four segments cho- 
sen to separate regions with and without strong NXB fea- 
tures. Overall, the spectra span the 0.7 — 10.5 (XIS03) 
and 0.7 — 7 (XIS1) keV energy ranges. Best-fit mod- 
els are found by minimizing the modification of the C- 
statistic [44[ implemented in Xspec. In addition to the 
astrophysical background components detailed above, 
additional NXB power-law and emission line components 
are independently included in the models in each segment 
- see [39| for details. The results of these fits are shown 
in Figure 2. 

Once these baseline models are established, we add an 
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FIG. 2. Unsubtracted Suzaku XIS03 {left panel, (a)) and XIS1 
(right panel, (b)) spectra of the Ursa Minor and Draco Dwarf 
Spheroidal Galaxies. The best joint-fit baseline models, that in- 
clude both non-X-ray and astrophysical background components, 
and residuals are also shown. The spectra are divided into seg- 
ments spanning the following energies: 0.7 — 1.3 keV (Ursa Mi- 
nor: black, Draco: light blue), 1.3 — 2.6 keV (Ursa Minor: red, 
Draco: magenta), 2.6 — 5.2 keV (Ursa Minor: green, Draco: yel- 
low), 5.2 - 10.5 keV for XIS03 or 5.2 - 7 for XIS1 (Ursa Minor: 
blue, Draco: orange). 



unresolved Gaussian emission line, stepped through the 
full bandpass in 25 eV steps, and derive 99% upper con- 
fidence levels (no significant detection was made) on the 
line surface brightness (that is permitted to be negative). 
The joint fits are conducted such that the three param- 
eters of interest, in addition to the line energy, are the 
sum of the Draco and Ursa Minor line surface bright- 
nesses and their ratio. Hence the 99% upper confidence 
level corresponds to A — statistic = 11.3. 



with an NFW profile with M 20 o = 3 - 30 x 10 9 Me@, 
and a scale radius determined by the WMAP5 mass- 
concentration relation [4?} . Within apertures of 7.7' 
- the region from which most of the source flux orig- 
inates - and for distances of 77 (Ursa Minor) and 76 
(Draco) kpc the corresponding dark matter surface 
density for a (tidal) truncation radius of 1.5 kpc [HI is 
~ 150 ± 50 M Q pc -2 for each galaxy, where the cutoff is 
implemented using the "n = 2 BMO" generalization of 
the NFW profile [49j . Including an additional summed 
Milky Way contribution of ~ 150 M Q pc -2 [H3|, we adopt 
a fiducial summed, total line-of-sight surface mass den- 
sity Sdm = 450 M Q pc~ 2 , mindful that there is an un- 
certainty of at least 100 M Q pc~ 2 . 




Dark matter surface densities 



0.2 0.4 0.6 
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FIG. 3. NFW mass profiles for M 20 o = 0.1, 0.3, 1.0, and 3.0 
xlO 10 Mq (solid, dotted, short-dashed, long-dashed curves, 
respectively) compared with estimates inferred by [i^] (filled 
symbols) and inferred or compiled by [ifjl ] (open symbols); 
triangles (squares) denote Ursa Minor (Draco). 



D. Constraints on moduli dark matter 

For our estimated summed dark matter surface mass 
density £dm = 450 M Q pc -2 , our 99% upper confidence 
levels on emission line surface brightness yields upper 
limits on the dark matter radiative decay rate over the 
Suzaku energy bandpass. We compare these limits as a 
function of energy (in 150 eV bins) for the joint XIS03 
analysis with the predicted decay rate of moduli dark 



Figure 3 shows the dynamically estimated mass profiles 
from [42j and [45[ , along with other determinations in the 
literature compiled by the latter. The Draco and Ursa 
Minor profiles are very similar within 600 pc, consistent 



where Af200 is the mass encompassing an overdensity, relative to 
the critical density, of 200 
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matter for b = 10 (equation 3) in Figure 4. Our up- 
per limits on b as a function of are shown in Fig- 
ure 5 for separate emission line limits from the XIS03 
and XIS1 joint spectral analysis. Significant constraints 
on the moduli dark matter candidate are obtained in 
the ^8 — 20 keV mass range from this initial inves- 
tigation. Significant improvements in sensitivity and 
mass range will be realized in observations made with 
the Astro-H Observatory^ scheduled for launch in 2014, 
with its capabilities for high energy resolution imaging 
spectroscopy in soft X-rays (Soft X-ray Spectrometer) 
and low background, moderate energy resolution, imag- 
ing spectroscopy in hard X-rays (Hard X-ray Imager). 



a line from decay of this relic dark-matter particles. If 



o 

0) 




2 4 6 

line energy (keV) 



FIG. 4. Comparison of excluded dark matter radiative decay 
rate (shaded histogram), and predictions for moduli decay 
assuming b = 10. The limits are based on 99% upper confi- 
dence levels on line fluxes from joint analysis of Ursa Minor 
and Draco Suzaku XIS03 spectra and a summed, line-of-sight 
dark matter surface mass density Edm = 450 Mq pc~ 2 . 



100 




<t> (keV) 



FIG. 5. Upper limit on the moduli b parameter represent- 
ing the (inverse of) the effective moduli scale in units of the 
reduced Planck scale, based on 99% upper confidence levels 
on line fluxes from joint analysis of Ursa Minor and Draco 
Suzaku XIS03 (red) and XIS1 (blue) spectra and a summed 
dark matter surface mass density Ed m = 450 Mq pc~ 2 . 



such a line is detected, the photon energy will provide 
information about the particle mass. 

We have presented the results of a search for such 
dark matter particles using the data from the first 
deep X-ray observations of the Draco and Ursa Minor 
dwarf spheroidal galaxies, dark matter dominated sys- 
tems with extreme mass-to-light ratios and low intrinsic 
backgrounds. The absence of an emission line results in 
new constraints on the relevant new physics. The limits 
thus obtained are sufficiently interesting so as to raise 
the prospect of a search with improved range and sen- 
sitivity - and possibly a detection - in the near future 
utilizing expected upgrades in the capabilities of X-ray 
observatories. 



VI. CONCLUSIONS 

String moduli and supersymmctry moduli produced in 
the early universe may persist as a form of dark matter. 
Current knowledge of supersymmetry breaking and the 
recent discovery of the 125 GeV Higgs boson, imply that 
the masses of moduli particles are above about a keV, 
and the most interesting range is between 1 keV and 100 
keV. In this range, X-ray telescopes may be able to detect 
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